Relationship between Entropic Bottleneck in Free Energy Landscape, Nonexponential 
Relaxation and Fragility of Glass-Forming Liquids 
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A mesoscopic model is proposed to explain the anomalous dynamics in a supercooled liquid as its 
glass transition temperature is approached from above. The model is based on the assumption of 
f3 organized a process, with the requirement of coherent excitation of a minimum critical number 
N c of j3 processes in the surroundings of a total Ng. Numerical evaluation of the model shows that 
the growth in this critical number in the background of a modest Np can lead to a severe entropic 
bottleneck and slow down the dynamics to an extent observed near real glass transition. The fragility 
of the glass-forming liquid is shown to be correlated with the growth of the ratio y(= N c /Ng). 



o 

•4-* 



I 

C 

o 
o 



> 

in 

m 
O 
m 
O 



T3 

O 
o 



PACS numbers: 64.70.Pf, 61.43.Fs, 66.20.+d, 



Relaxation in a viscous liquid above but near its glass 
transition temperature T„ has been a subject of immense 
attention in recent time 0, Li Lj [J bl L* • The dramatic 
increase in the characteristic relaxation time r (which 
varies as much as 15 orders of magnitude for fragile liq- 
uids 0, @) is most frequently described by the Vogel- 
Fulcher-Tammann (VFT) equation [t| 



t = A vft exp[B/(T-T )], 



(1) 



where Ay > B and T are found to remain constant over 
a range of 4 — 5 orders of magnitude variation of t. To, 
known as the Vogel temperature, is typically 30—50 K be- 
low T g . Tq = corresponds to the Arrhenious behavior, 
characterizing the strong liquid limit. The fragile behav- 
ior is further marked by nonexponentiality in the relax- 
ation functions, which can be described by a stretched 
exponential or Kohlrausch- Williams- Watts function 



0(t) = exp[-(t/T) 



01 



(2) 



where /3, < < 1, is the stretching exponent. For a 
typical fragile glass-former, [3 generally decreases from 
near 1 at high temperatures to below 0.5 close to T g 
with a display of nearly monotonic temperature depen- 
dence Experimental and computational studies sug- 
gest that the stretched exponential relaxation may owe 
its origin to the growth of spatially heterogeneous re- 
laxing domains, where each individual exhibits a nearly 
exponential relaxation with a relaxation time that varies 
significantly among the individuals 0, . Such hetero- 
geneous domains have been found experimentally to span 
2 — 3 nm [l2l[T3| . However, the origin of this modest size 
of the heterogeneous domains is not clearly understood at 
present. A consistent description of all the above aspects 
of dynamics in the supercooled regime has remained a 
major scientific challenge over many decades. 

Several theories suggest the appearance of a growing 
length scale as the glass transition is approached from 
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above [TJ|. The celebrated Adam-Gibbs (AG) theory 
|15| attempts to explain the temperature dependence 
of relaxation time in terms of a temperature variation 
of the size z of the cooperatively rearranging region 
(CRR); the lower limit z* is "shown" to follow the re- 
lation z* = Nas*/S c , where s* is the critical configura- 
tional entropy of a CRR corresponding to z*, S c is the 
molar configurational entropy of the macroscopic sam- 
ple, and Na is the Avogadro constant. According to AG 
theory, the sluggishness near T g is due to the scarcity of 
the number of configurations available to a CRR, which 
is reflected in the increasingly smaller value of S c ; the 
resultant rapid increase in the value of z* leads to the di- 
vergence of the relaxation time r. However, experimental 
[l2| and computer simulation 0] studies have failed to 
find a convincing evidence of a growing length scale near 
T g . Nevertheless, AG theory of entropy crisis provides a 
useful conceptual framework. 

Perhaps the most successful quantitative theory of re- 
laxation phenomena in the supercooled liquid state is the 
mode-coupling theory (MCT) 0, . It is known to be- 
come inadequate at low temperatures below T c , called the 
MCT critical temperature. This is presumably because of 
the prevalence of the thermally activated hopping 
below T c (unaccounted for in MCT) between the adjacent 
minima of the energy landscape (i. 2QJ . Recent computer 
simulation studies have revealed that hopping is a highly 
cooperative phenomenon, promoted by many body fluc- 
tuations |2l|. Computer simulations have further shown 
that large amplitude hopping of a tagged particle is of- 
ten preceded by somewhat larger than normal, still small 
amplitude, but collective motion, of its neighbors. A sim- 
ilar picture of (3 and a processes was earlier proposed by 
Stillinger, where the (3 process corresponds to transitions 
between inherent structures (IS) within the meta-basin 
while a process corresponds to transitions between deep 
meta-basins 6]. Motivated by these findings, we here 
present a mesoscopic model of relaxation in supercooled 
liquids (~ below T c ), where an a process is promoted by 
coherent excitations of a minimum number of (3 processes 
within a CRR. The (3 processes are assumed, as a first 
approximation, to occur independently. The requirement 
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of coherence is that a given minimum number N c among 
the total number Np of (3 processes must be in the ex- 
cited state during a small interval for an a process (i.e., 
transition out of a meta-basin) to occur. In our model, 
a CRR is characterized by an Np number of identical 
noninteracting two- level systems (TLSs). Each of these 
TLSs transits back and forth at equilibrium at tempera- 
ture T between its two levels, labeled (ground) and 1 
(excited). The two levels are separated by an energy e. 
The waiting time before a transition can occur from the 
level i(= 0, 1) is random, but is drawn from a Poissonian 
probability density function given by 



tpi(t) = —expi—t/n), 
n 



0,1, 



(3) 



(4) 



where r, is the average time of stay in the level i. If pi 
denotes the canonical equilibrium probability of the level 
i being occupied, detailed balance gives the following re- 
lation 

Po t 

where K is the equilibrium constant for the two levels. 

We define a variable C, 3 (t), (j = 1, 2, , Np), which takes 

on a value if at the given instant of time t the level 
of the two-level system j is occupied and 1 if otherwise. 
Cj(t) is thus an occupation variable. The control variable 
Q(t) is defined as 



(5) 



3=1 



Q(t), which serves as an order parameter, is a stochastic 
variable in the discrete integer space [0, Np] and carries 
information of the excitation prevailing in the CRR at 
time t. The rate of a relaxation depends crucially on Q. 
For simplicity, we assume here that the relaxation occurs 
with unit probability at the instant Q reaches N c , an 
integer greater than the most probable value of Q, for 
the first time. This restriction can be removed, but only 
at the expense of the analytical solution presented below. 

We solve the model for the average relaxation time r, 
for a a given pair of N c and Np, by using the method 
of mean first passage time j23| • The probability that the 
stochastic variable Q takes on a value I at time t, P(l; t), 
satisfies the following master equation 



dP{l;t) 
Jt 



[(Np-l + l)/T ]P(l-l;t) 

+ [(l+l)/ n ]P(l + l;t) 
- [(Np-l)/r ]P(l;t)-(l/n)P(l;t). (6) 

The mean first passage time t(V), which is the mean time 
elapsed before the stochastic variable Q (starting from 
its initial value Z < A^ c — 1) reaches N c for the first time, 
satisfies the following equation related to the backward 
master equation: 

[(iV^-0/To][r(Z+l)-r(0] + a/n)[ra-l)-ra)] = -1, 

(7) 
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FIG. 1: Semi- log plot of the single domain relaxation function 
Sso(t) versus t at a given T. The figure is based on the data 
obtained by simulating the model for single domain relaxation 
with e = ksT, \ib being the Boltzmann constant, and Np = 
10 for N c — 6, 7 and 8 (dotted line, long dashed line, and 
dot-dashed line, respectively) . The solid lines are the straight 
line fits. 



subject to an absorbing boundary condition at Q — N c , 
t(N c ) = 0, and a reflecting boundary condition at Q = 0, 
t(— 1) = r(0). We solve Eq. (7) to obtain t(1) as a sum 
over hypergeometric functions F(a, b; c; z) 



-(I) = T (l + 1/K) 



N/3 



N c -1 

E 



FjNp + 1, Np - n; Np - n + 1; -l/K) 
Np-n 



(.8) 



Eq. (8) shows that when the two states are of the same 
energy, that is, K — 1, even then the relaxation slows 
down significantly. This is purely an entropic effect - a 
nice example of entropic bottleneck [2^ . When there is 
an energy bias against the excited state (state 1), the 
bottleneck becomes more severe, as shown below. 

The relaxation of the relevant time correlation func- 
tions, such as stress and density, is assumed to be caused 
by the a relaxation. The time dependence of the a 
relaxation can be quantified by the survival probabil- 
ity correlation function Ssd, of the initial meta-basin 
state. This is obtained by first simulating the model for 
single domain relaxation. The parameter 7, defined as 
7 = N c /Np, provides a measure of the free energy bar- 
rier to relaxation having both the energy and the entropy 
contributions. Fig. (1) plots the time dependence of the 
single domain relaxation function Ss d {t) in the logarith- 
mic scale at a given temperature for three values of N c . 
The straight line fits show a single exponential decay of 
Sso(t)- The time constants obtained are identical to 
those obtained from our analytical expression, given by 
Eq. (8). The single domain relaxation is found to slow 
down considerably, as expected, with increasing 7. Fig. 
(2) illustrates an exponential dependence of the single 
domain average relaxation time t$d on Np at fixed 7 
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FIG. 2: Dependence of the single domain average relaxation 
time tsd on Np at a given temperature T with e/(fcgT) = 1 
for different fixed values of 7. The inset shows the 7 depen- 
dence of the slope m of the semi- log plot of the scaled tsd 
versus Ng. 
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FIG. 3: The long time behavior of the averaged relaxation 
function S(t). The solid line is a fit to the stretched exponen- 
tial function with r = 1.100 x 10 6 n and j3 = 0.670. The inset 
shows the temperature dependence of the stretching exponent 
f3. 



and T; t sd is a weighted average of r(Z) over the initial 

N c -1 

distribution: tsd = ^ P(1;0)t(1). The dependence, 

1=0 

characterized by the slope of the semi-log plot of the 
scaled tsd versus N/3, becomes stronger as 7 increases. 
The inset suggests a power law behavior of the slope m, 
m ~ "f v with v ~ 3.4. Note that tsd is also the mean 
waiting time. 

In a heterogeneous environment within a bulk sam- 




FIG. 4: Arrhenious plot showing the scaled characteristic re- 
laxation time t as a function of reduced inverse temperature 
Tg/T. The two sets of data, marked with circles and squares, 
for r are obtained from the stretched exponential fit of the 
long time behavior of the average relaxation function S(t). 
The circles correspond to the data set obtained from a linear 
variation of 7 from 0.6 to 1.0 as T g /T is increased from 0.68 to 
1.0 with N/3 held fixed at 20. The squares denote the data set 
corresponding to a linear growth of 7 with T g /T, both rising 
from 0.70 to 1.0, with iV c and Ng varied together, the latter 
from 10 to 20. The solid and the dotted lines, which are the 
corresponding fits to the VFT form with To = 0.689 T g and 
0.789 T g , respectively, illustrate fragile liquid behavior. The 
data set, marked with triangles, are obtained with 7, Ng, and 
e all held constant: 7 = 1.0, Ng = 40, and e = fcsT m . The 
data are scaled so as to have the same r value at T g as for 
the other two sets. With no distribution of e, the relaxation 
is exponential with the time constant r. The straight line 
fit corresponding to the Arrhenious behavior is typical of a 
strong liquid. Here, for the sake of simplicity, numerical eval- 
uation is done by taking tsd = r(Q rap ), where Q mv is the 
most probable value of Q; this involves negligible error. 



pie, a fluid-like region can be characterized by having, 
at a given time, on the average, a relatively large num- 
ber of (i processes in the excited state. The reverse is 
true for solid-like regions. Therefore, a Gaussian distri- 
bution of e among CRRs can incorporate the existence 
of heterogeneous domains into the model (one can also 
assume an exponential distribution). The present calcu- 
lation takes the mean < e > of the distribution to be 
unity and the standard deviation a = 0.05 in the units 
of ksTm, T m being the melting temperature. T m is also 
used to define T g : T g = 2T m /3. Furthermore, the model 
assumes 7V C to grow as the reduced inverse temperature 
Tg/T increases until N c reaches Np at T g . The Gaussian 
distribution of e results in a continuous distribution of 

tsd- The average relaxation function S(t) is calculated 

00 

from S(t) = J dT S Dg{TsD)exp(-t/T S D), where g(T S D) is 


the probability density function. The long time behavior 
of S(t) fits well to the stretched exponential function as 
shown in Fig. (3). The inset shows a monotonic decrease 
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of the stretching exponent (3 with T g approached from 
above, as indeed observed experimentally. The study of 
the temperature dependence of the characteristic relax- 
ation time t requires time to be scaled by the same unit 
at all temperatures. t± at T g , which we denote by ri )S , is 
chosen for that. The scaling needs, for a transition state 
theory (TST) calculation, an input to , the energy bar- 
rier to the transition from the level to 1, which is taken 
to be 4 ksTm and held constant. Fig. (4) shows the 
T g scaled Arrhenious plot of the characteristic relaxation 
time. The fit to the VFT equation is good. The growth 
of 7 with T approaching T g from above determines the 
fragility, as illustrated in Fig. (4). Strong liquid limit 
is obtained with 7 held fixed at 1. Strong liquids are 
spatially (and, dynamically) correlated even at high tem- 
peratures. This justifies the choice of a higher 7 and Np 
values in the strong liquid limit. 

We have also calculated the waiting time distribution 
W(t) (WTD) |2J| for this a process. Even when the 
distribution of energy gap in the meta basin is Gaussian, 
WTD is non-Gaussian with a stretching at long r. The 
waiting time distribution gets modified if the distribution 
of energy gap (e) is exponential - the most affected region 
is obviously the small r limit. 

The present model can be considered as a representa- 
tive in a simple form of a class of wider, more general 
models. An immediate generalization will be to include 



the correlations among the (3 processes within a CRR. 
Note that the present model explains the slow down of 
relaxation in fragile liquids as result of an entropic bot- 
tleneck superimposed on the energy constraint and does 
not require any diverging length scale. The model pre- 
dicts that the observed 13 orders of magnitude increase 
in relaxation time originates from the combined effect of 
energy and entropy, and it is not possible to separate the 
two effects. The model is based on a requirement of dy- 
namical correlation between the (3 and a processes. This 
correlation itself depends on the time scale of separation 
between the two. In the language of energy landscape 
@ IH. E3 > the (3 transitions assumed here occur within 
the super-structures (meta-basins) of a deep minimum 
and they promote transitions between two deep minima - 
that is, among the meta basins. While the model rests on 
dwindling entropy as the temperature is lowered, it does 
not invoke any thermodynamic phase transition. Most 
notably, only a modest growth in the size of CRR, repre- 
sented here by the size of Np, is required to capture the 
experimentally observed slow down. 
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